Effect of unitary impurities in non-STM-types of tunneling in high- Tc superconductors 
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Based on an extended Hubbard model, we present calculations of both the local {i.e., single- 
site) and spatially-averaged diflterential tunneling conductance in d-wave superconductors containing 
nonmagnetic impurities in the unitary limit. Our results show that a random distribution of unitary 
impurities of any concentration can at most give rise to a finite zero-bias conductance (with no 
peak there) in spatially-averaged non-STM type of tunneling. This is in spite of the fact that local 
tunneling in the immediate vicinity of an isolated impurity does show a conductance peak at zero 
bias. We also find that to give rise to even a small zero-bias conductance peak in the spatially- 
averaged type of tunneling the impurities must form dimers, trimers, etc. along the [110] directions. 
In addition, we find that the most-recently-observed novel pattern of the tunneling conductance 
around a single impurity by Pan et al. [Nature 403, 746 (2000)] can be explained in terms of a 
realistic model of the tunneling configuration which gives rise to the experimental results reported 
there. The key feature in this model is the blocking effect of the BiO and SrO layers which exist 
between the tunneling tip and the Cu02 layer being probed. 



PACS numbers: 74.50. -f-r, 74.62. Dh, 74.80.-g 
I. INTRODUCTION 

Several years ago, one of us (CRH)0 showed that 
the quasi-particle spectrum of a c?-wave superconductor 
(DWSC) contains a special class of excitations — called 
midgap states (MS's) in that reference — which have es- 
sentially zero energy with respect to the Fermi energy, 
and are bound states with their wave functions localized 
at the vicinities of various kinds .-pf, defects rin the sys- 
tem, such as a surfaceEJ and a grainaa or twiuEl boundary. 
These MS's form an "essentially dispersionless" branch of 
elementary excitations, in the sense that their momenta 
along a flat surface or interface can essentially range from 
—kp to kp (Fermi momentum), and yet with almost no 
accompanied kinetic energy variation. The existence of 
surface and interface MS's appears to have^ already been 
confirmed by several types of experiments.ETu An impor- 
tant question is whether a unitary impurity can also give 
rise to MS's. The MS's have a topological origin, in the 
sense that in the semi-classical WKBJ approximation, 
which makes these states truly midgap, their existence 
requires the satisfaction of one or more sign conditions 
only. More precisely, one can describe such a state in 
terms of one or more (as a linear combination) closed clas- 
sical orbits, each of which must encounter two Andreev 
reflections by the pair potential at two different points 
of the Fermi surface where the pair potential (A(/c)) has 
opposite signs. For each MS formed at a specular surface, 
only one such closed i^rbit is involved, so only one sign 
condition is required.EJ For a MS formed at a flat grain 
(or twin) boundary, modeled as a planar interface with 
transmissivity < t < 1 and different crystal orientations 
on its two sides, two closed classical orbits are involved, 
corresponding to the possibility of transmission and re- 
flection at the interface. Thus two sign conditions must 



be satisfied simultaneously.§ If a unitary impurity could 
be represented by a circular hole of a radius much larger 
than the Fermi wavelength, then quasi-classical argument 
can be expected to hold, and at least some non- vanishing 
number of midgap states should exist near its boundary. 
But when a unitary impurity is of atomic size, such a 
quasiclassical argument becomes dubious. In fact, if one 
thinks of scattering by an impurity as a linear combina- 
tion of infinite number of classical orbits, corresponding 
to the possibility of scattering by all angles on the Fermi 
surface, then it would seem that an infinite number of 
sign conditions would have to be satisfied in order for 
the impurity to induce some midgap states, correspond- 
ing to requiring sign change of the pair potential between 
any two points on the Fermi surface, which clearly is not 
satisfied. This argument suggests that midgap states can 
not form near an impurity, or at least they will not be 
exactly midgap even in the semiclasscal approximation. 
However, this argument does not distinguish between a 
unitary impurity and a non-unitary one. On the other 
hand, early treatments of a random distribution of uni- 
tary impurities in DWSC's, based on the self-consistent 
t-matrix approximation, have indicated a broajd peak 
at zero energy in the density of states (DOS).lj (Even 
earlierj-aimilar studies on p-wave SC's also showed such 
peaksJlj) Perhaps Balatsky et al. are the first to men- 
tion within this context that a single unitary impurity in 
a c?-wave SC will lead to twO|-aero-energy bound states 
(per spin), with energies ieo.til (Buchholtz and Zwick- 
nagl made a siipilar statement much earlier for p-wave su- 
perconductors J13) The t-niatrix approximation employed 
in all of these works, which includes a semiclassical ap- 
proximation, restores particle-hole symmetry, which is 
not a very good approximation in high- Tc superconduc- 
tors (HTSC's). Without this symmetry the energies of 
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these "MS's" are most-likely not exactly zero. Unlike 
MS's at surfaces and interfaces, which involve a small 
number of points on the Fermi surface only, and there- 
fore can easily avoid the gap-node directions, in the case 
of a state localized around an impurity, all points on the 
Fermi surface must be involved, including the nodal di- 
rections where the gap vanishes. Then for a state whose 
energy is not exactly zero, its wave function must be 
able to leak to infinity near the nodal directions. That 
is, such a state can not be a genuine bound state, and can 
only be a resonant state. When there is a finite density 
of impurities, such wave functions should then possess 
a long-range interaction with each other via the "leak- 
ages" , leading presumably to a broad "impurity band" . 
The self-consistent f-matrix approximation does not take 
this point into account properly^^nd is therefore not sat- 
isfactory. (However, see Joynlllj for an opposite view, 
except on the assumption of particle-hole symmetry and 
the validity of neglecting "crossing diagrams" when per- 
forming impurity averaging in two dimensions. But we 
think that neglecting the crossing diagrams is precisely 
why that approximation can not treat the impurity in- 
teraction properly.) A numerical|-approach has been in- 
troduced by Xiang and WheatleyEJ to avoid these short- 
comings. For a single unitary impurity and a random 
distribution of unitary impurities it gave results in good 
agreement with the t-matrix approximation on the DOS. 
Numerical solutions of a lattice BCS model with nearest 
neighbor attraction has also been employed by Onishi et 
a/.Jlj who showed that: (i) their exist two essentially- 
zero-energy states (per spin) localized around an impu- 
rity, (ii) their wave functions have long tails in the nodal 
directions, (iii) as a result such states localized around 
two impurities separated by a large distance in compar- 
ison with the coherence length can still interact with 
each other, leading to a broad impurity band. Thus the 
general picture outlined above appears to be confirmed. 
But Ref. |lj probably also assumed particle-hole symme- 
try since it plotted the DOS for positive energy only, 
whereas in Ref. |l5| only ^ = is considered which has 
exact particle-hole symmetry. Numerical diagonaliza- 
tion method has also been applied to particle-hole non- 
symmetric models to study impurity effects, but without 
addressing the questions raised here. See, for example, 
Ref. Rather, that work and several other worksEZi 
debated on whether there is localization in the impurity 
band — a topic which is not our concern here. Instead, 
we wish to address the roles played by unitary impu- 
rities in single-particle tunneling. In particular, it has 
been noted that MS's formed at surfaces and interfaces 
of DWSC's can lead to an obserpjahle psro-bias conduc- 
tance peak (ZBCP) in tunnelingBclEa til Indeed, several 
recent tunneling experiments pierformed withEj STM/S 
and other tunneling techniquef£j on HTSC single crys- 
tals and epitaxial thin films have shown that in a6-plane 
tunneling a very prominent ZBCP can be observed, es- 
pecially on {110} surfaces. The observed ZBCP exists 
continuously for long distances along the surfaces, and 



the observed tunneling characteristics can be quantita- 
tively fitted by a generalized Blonder-Tinkham-Klapwijk 
theorynJ which includes the effects of the MS's formed 
on such surfaces of DWSC's. Assuming that this in- 
terpretation is correct, a question one can ask next is: 
Can unitary impurities be responsible for at least some of 
the observed ZBCP's? This question is meaningful since 
ZBCP's have been observed ubiquitously in all kinds of 
tunneling settings, and some of them may not possess 
surface and interface MS's. (See Ref. ^ for a review.) 
Theoretical results reviewed above, although not conclu- 
sive, seem to suggest that unitary impurities can also 
give rise to ZBCP's. But experimental evidence in non- 
STM/S types of tunneling seem to suggest the contrary. 
To see that this is the case, one must first exclude ZBCP's 
observed in ab-plane tunneling, and on tunneling per- 
formed on polycrystaline and ceramic samples, since in 
these cases contributions from the surface and interface 
MS's can most-likcly dominate. (Whereas no MS's can 
form on a flat {100} surface of a_pWSC,lil surface rough- 
ness can reverse this conclusion.Ej Thus we exclude a-axis 
and {nOm}-directional tunneling as well, if n=f 0.) One is 
left with c-axis tunneling on single crystals and epitaxial 
thin films only. (Nominal c-axis point-contact measure- 
ments may actually be seeing some a6-plane tunneling, 
since the tunneling tips have been pushed into the HTSC 
samples in these measurements. So they should also be 
excluded.) Even in the surest c-axis tunneling cases one 
must still distinguish between STM/S- type tunneling, 
which can explore the tunneling characteristics in the 
close vicinity of a single isolated impurity, and the other 
non-localized tunneling techniques, which see a spatially 
averaged tunneling characteristics. For the later type, 
the spectral weight of the impurity contribution must not 
be too low to be observed, so the impurity concentration 
must not be too low. It is this kind of tunneling which 
we are interested in here, since we suspect that once the 
spectral weight is sufficiently large, the interaction be- 
tween a random distribution of impurities will also be so 
large that it still can not give a ZBCP, but only a finite 
conductance at zero bias, G(0), as a local minimum or 
even an extra dip. Searching the literature, we find three 
more recently published papers reporting the STM/S re- 
sults on the observatiatL q£ a ZBCP-like feature in the 
vicinity of an impurity.tjLj For non-localized tunneling, 
we find at most a few cases which can weakly suggest 
that the small ZBCP's observed in them might originate 
from impurities.Ej (Even in these cases, it is not clear 
whether MS's could have formed at some exposed Cu02 
edges at the interface with the insulating barrier. The 
epitaxial films might also have grain boundaries which 
could host MS's.) Most c-axis tunneling data which ex- 
hibit a clear gap feature show a minimum at zero bias, 
with some showing essentially simple d-wave behavior, 
with very small G(0),l3 and sorne^howing nearly d-wave 
behavior but with a finite G(0)EJ Still others show fea- 
tures on both sides of zero bias, givipg the impression 
that there is an extra dip at zero bias.cHi Very few publi- 
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cations seem to have systematically studied the impuritjj 
effects in c-axis tunneling. We find one: Hancotte et alSn 
showed that Zn substitution (< 1%) in the Cu02 planes 
of BSCCO (2212) caused G{0) to markedly increase, ac- 
companied by a reduced gap. Here the spectral weight 
of impurity effects is clearly large enough, yet not even a 
trace of a ZBCP was observed. 

More recently, the quasiparticle properties around a 
single impurity have been imzKaiigated in more detail in a 
two dimensional t-J model,LjO with a focus on whether 
the ZBCP observed with STM/S near an impurity is split 
or not. The purpose of this paper is, on the other hand, 
to attempt to answer the following precise question: As- 
suming that HTSC's are DWSC's, and isolated unitary 
impurities do possess near-zero energy resonant states 
which can be observed as a ZBCP-like feature by STM/S 
in the close vicinity of such an impurity, can any concen- 
tration of a random distribution of unitary impurities be 
able to give rise to an observable ZBCP (of any width) 
in non-STM/S types of tunneling, or there must be some 
spatial correlation in the impurity distribution before a 
ZBCP can appear in such types of tunneling? For this 
purpose we have performed an extensive numerical study. 
We have introduced the supercell technique so that the 
finite-size effects from the exact diagonalization can be 
overcome, and the desired energy resolution can be ob- 
tained. This technique has the ability to treat well the 
impurity of atomic size. Moreover, the band structure 
effects can be incorporated in a natural way. The results 
show the conductance behavior is sensitive to the posi- 
tion of the chemical potential within the band and to the 
impurity configuration at the atomic scale. Our results 
have indeed confirmed our suspicion that for a simple 
random distribution of unitary impurities, either their 
spectral weight is too low for their effects to be observ- 
able in non-STM types of tunneling, or their interaction 
is so strong that only a finite G{0) is obtained as a local 
minimum rather than a peak, because the impurity band 
has spread wide, with the center of its contribution to the 
density-of-states function lower than its two sides. In ad- 
dition, we find that if only enough number of the unitary 
impurities form nearest neighbors along the [11] direc- 
tions in a Cu02 plane, henceforth called "[ll]-directional 
dimers", then a weak ZBCP can appear in non-STM 
types of tunneling. (We also find that if enough such 
impurities form [ll]-directional trimers, then the ZBCP 
can be even taller and narrower, but the chance of form- 
ing such alignments in an actual sample is probably very 
low. On the other hand, we think that dimers can proba- 
bly form with not very low probability.) Furthermore, we 
have also shown that the most recently observed pattern 
(by STM/S) ofpthe local tunneling conductance around a 
single impurityEa can be explained by taking into account 
that the STM/S tip in that experiment is separated from 
the Cu02 plane under probe by a BiO layer and a SrO 
layer. Therefore, the tunneling tip can not communicate 
with the atom directly below it in the Cu02 plane, due 
to the blocking effect of the atoms directly above it in the 



BiO and SrO layers. Rather, we think that the measured 
"local tunneling conductance" by the tunneling tip on a 
Cu site in the Cu02 plane is actually that averaged over 
a local region around that site, excluding the contribution 
from the atom at that site, because of the blocking effect 
just described. Using this very reasonable postulate, we 
find that we can at least semi-quantitatively understand 
the pattern observed in Ref. |2^, including why it peaks at 
the impurity site, and vanishes at its four nearest neigh- 
bor Cu sites, etc. 

II. THEORETICAL METHOD 

To model decoupled copper-oxygen layers in HTSC's, 
we consider the single-band extended Hubbard model de- 
fined on a two-dimensional square lattice (lattice con- 
stant a) with nearest-neighbor hopping, and on-site vf^ 
pulsivc and nearest-neighbor attractive interactions .Eil 
For our purpose, we introduce supercells each with size 
N^a X Nya. We then define the supercell Bloch states 
labeled by a wave vector k and a site index i within the 
supercell. In the mean field theory, the task becomes 
to exactlyrdiagonalize the Bogoliubov-de Gennes (BdG) 
equations :E3 

_/iJy(k) Ay(k) \ / Ul'-X (ul^\ 

At.(k) -ijy(k) ; 1^ ) - ^-'^ 1, .--^ ) ■ 

(1) 

Here u^'^ and w"'*' are the Bogoliubov amplitudes corre- 
sponding to the eigenvalue i?„,k- 

iJij (k) = -te^'-^'^6,^s.j + (Ui - t^)Su , (2) 

Ay = Ao(i)^ij+A^(i)e^'^ '^'^5;^^ . , (3) 

where t is the hopping integral, fi is the chemical po- 
tential, 6 = ±x, ±y are the unit vectors along the crys- 
talline a and b 2Trnx.y/Mx,yNx,ya with 
nx,y — 0, 1, 2, . . . , M^^^y — 1. Mx^yNx^ya is the linear 
dimension of the whole system, which is assumed to 
be made of M = x My super-cells. The single- 
site nonmagnetic impurity scattering is represented by 
Ui — CA) X^i'e/ with the summation over the set of 
impurity sites. The self-consistent pair potentials are in 
turn expressed in terms of the wavefunctions {ui^Vi): 

^o(i) ^ 17 E "r^^^'"^)* tanh(i?„,k/2fcBr) , (4) 

n.k 

ri,k 

X tanh(S„,k/2fcBT) , (5) 
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where fc^ is the Boltzmann constant and T is the abso- 
lute temperature. In the case of repulsive on-site (50 < 0) 
and nearest-neighbor attractive [gi > 0) interactions, the 
d-wave pairing state is favored, the amplitude of which 
is defined as: Arf(i) = [Axji) + A_i„(i) - Ax,(i) - 
A-x,(i)]/4. 

We define the local differential-tunneling-conductance 
(DTC) byEll 

+ K''fnE^,u + E)], (6) 

where the prefactor 2 comes from the two-fold spin de- 
generacy and f'{E) is the derivative of the Fermi distri- 
bution function f{E) = [expiE/ksT) + Gi{E) can 
theoretially be measured by STM experiments. (How- 
ever, see later for a possible complication when the 
probed layer and the STM/S tip are separated by other 
atomic layers.) The spatially-averaged DTC is defined 
by G{E) EE Y>iGi{E)/NxNy, which can essentially be 
measured in many non-STM/S types of tunneling exper- 
iments (planar, ramp, etc.). (But squeezable junctions 
might measure something in between, depending on the 
pressure applied.) 

III. NUMERICAL RESULTS 

In the numerical calculation, we take the supercell size 
Nx = Ny = 35, the number of supercells M = 6 x 6, and 
the temperature fc^T — 0.02^. The lattice sites within 
one supercell are indexed as {ix,iy) with and iy each 
ranging from 1 to 35. In addition, we use the bulk value 
of the order parameter as an input to diagonalize Eq. (|l|) 
and the deformation of the order parameter near the im- 
purity is ignored. (This approximation should be quite 
acceptable for the questions we wish to answer here.) In 
the bulk system, the d-wave order parameter has the form 
Ak = 2Arf( cosk^a — coskyO), where k^^y are the x- and 
y-components of the wave vector defined on the whole 
system. Table | lists the resulting d-wave pair potential 
Ad and coherence length ^0 = hvp/T^Amax for the cho- 
sen values of 171 and /i, where is the Fermi velocity and 
Amax = 4 Ad is the maximum energy gap. The choice of 
the maximum energy gap is consistent with the bulk gap 
structure exhibited in the bulk DOS. In the calculation 
of Ad, an arbitrary negative value of go can be taken. 

A. The cases of a single impurity and a small cluster 
of impurities 

The quasiparticle property near a single impurity in a 
DWSC is complicated, and whether the energy of quasi- 
particle resonant states is exactly zero (relative to the 
Fermi energy) or whether the zero (or near-zero) energy 



states are split is very sensitive to the band structure and 
the impurity strength. To give a clear answer or clue to 
this question, we first consider the case of weak or mod- 
erately strong impurity. Figure |^ plots the local DTC as 
a function of bias directly on the single-site impurity (a) 
and one lattice constant away (b). As shown in Fig. [l|, 
when the impurity scattering is weak (Uq = 2.5t), the 
local DTC on the impurity site has a peak below the 
Fermi energy, which is consistent with the earlier study 
within the continuum theoryjlil while that at the site 
nearest-neighboring to the impurity has a double peak 
structure, one above and the other below the Fermi en- 
ergy. When the impurity scattering becomes stronger 
[Uq — lOi), the peak on the local DTC on the impurity 
site is pushed toward the Fermi energy with the ampli- 
tude strongly suppressed; nevertheless the double peaks 
on the local DTC at the nearest-neighbor site converge to 
each other and the intensity is enhanced. Since /i = in 
this calculation, the band is globally particle-hole sym- 
metric. Therefore the non-zero energy resonant states 
shown above originates from the local particle-hole sym- 
metry breaking. Differently, when the impurity scatter- 
ing goes to the unitary limit, whether there exists the 
particle-hole symmetry depends on solely on the position 
of chemical potential within the band. 

In all the following discussions, the single-site impu- 
rity potential strength is taken to be Uq = lOOi, so it 
is practically in the unitary limit. In Fig. || the local 
DTC is plotted as a function of bias near one single- 
site impurity located at site (18, 18) for various values 
of fj, but with gi — 2t fixed. The bias is normalized to 
Amax- The tunneling point is one lattice constant away 
from the impurity along the (10) direction, i.e., at (19, 
18). As is shown, when = where the particle- hole 
symmetry holds, a sharp zero-bias conductance peak is 
exhibited. When fi deviates from zero, the line shape 
of the local DTC becomes asymmetric with respect to 
the zero energy position. The peak near zero energy 
is seen to show an incomplete splitting and the height 
of this peak is decreased. The extent of splitting in- 
creases with the deviation of fi from zero. This result 
is quite different from the case of a {110}-oriented sur- 
face, where an un-split zero-energy peak in the local DTC 
shows up regardless of the position of the chemical po- 
tential. The similarity between the result obtained here 
by solving the extended Hubbard model and the corre- 
sponding result we obtained earlier by solving the t-J 
modclEj indicates that our conclusion about the splitting 
of the ZBCP is, at least in the mean field level, model 
independent. Next, we show that when more than one 
unitary impurities are present in the SC, the local DTC 
heavily depends on the impurity configuration. In Fig. 
we plot the local DTC in a DWSC with two and three im- 
purities forming a [ll]-directional nearest-neighbor dimer 
and trimer, respectively, for gi = 2t and fi — —t. With- 
out loss of generality, the impurity positions are placed 
at sites (17,19) and (18,18) for the dimer case, and at 
sites (17,19), (18,18) and (19,17) for the trimer case. The 
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measured position is still at site (19,18). As can be seen 
clearly, the near-zero-energy peaks in the local DTC for 
the single impurity case are strongly pulled toward zero 
energy in the dimer case. In the trimer case, the zero- 
energy peak is even more pronounced. To check this 
point more seriously, we have also calculated the low- 
est positive eigenvalues for the system composed of one 
supercell with one single impurity, and two impurities 
with different relative positions for gi = 2t and fi = —t. 
We find that a single impurity leads to two near-zero- 
energy states per spin, one just above and one just below 
the Fermi energy. The absolute energy corresponding to 
these two eigenstates is roughly Q.lt. (It is not zero be- 
cause there is no particle-hole symmetry in this model 
study.) But for the two impurities positioned as nearest 
neighbors of each other along the [11] direction, the low- 
est positive eigen-energy is only about O.OSi. If two im- 
purities are positioned as nearest neighbor of each other 
along the [10] direction, the lowest positive eigenvalue is 
roughly 0.09<, which is very close to the value for isolated 
impurities. In Fig. ^, we have plotted the variation of the 
lowest positive eigenvalue with the distance between two 
imputities aligned along the [11] direction. When two im- 
purities are far apart, the lowest positive eigenvalue os- 
cillates between ~ O.lOSt and '-^ 0.079t, around the value 
^ Q.lt for isolated impurities, as the distance between the 
two impurities are increased. This oscillation is a clear 
indication of the long range interaction between two im- 
purities that we have already discussed. (The long tails 
of the wave function of a bound state around a unitary 
impurity in the nodal directions are oscillating in Fermi 
wavelength, so the sign of the interaction between two 
impurities can change with distance.) When two impu- 
rities form a [ll]-directional dimer, they combine to play 
the role of a very short {11} edge. That is, they drasti- 
cally enhanced the probability for specular reflection by 
this edge, so the nearest-to-zero eigen-energies become 
very close to zero energy. This analysis demonstrates 
that the closest-to-zero eigen-energies due to the pres- 
ence of impurities are very sensitive to the short-range 
correlations in the impurity configuration. Furthermore, 
we plot in Fig. ^ the number of near-zero eigen-energy 
states versus the number of impurities aligned consecu- 
tively along the [11] direction. In this calculation, each 
state with its eigenvalue smaller than 0.04t is counted as 
a near-zero-energy state. As the result shows, the num- 
ber of near-zero-energy states increases almost linearly 
with the length of a [ll]-directional impurity line. We 
also find that, as the number of impurities is increased, 
there are more and more states with their eigen-energy 
approaching zero, which confirms from an alternative as- 
pect that the appearance of the ZBCP in the {110}- 
oriented DWSC tunnel junctions is independent of the 
position of ^ in the band. ■— ■ 

Most recently. Pan et aZr3 have performed an STM 
study on the effect of a single zinc (Zn) impurity atom 
on the quasiparticle local density of states (LDOS) in 
BSCCO. Besides revealing the predicted highly local- 



ized four-fold quasiparticle cloud around the impurity, 
the imaging also exhibited a novel distribution of the 
near-zero-resonant-energy LDOS near the impurity. In 
contrast to the existing theories, which give a vanishing 
LDOS directly at the site of the unitary impurity, and 
vanishing or low values on all atomic sites along the [11] 
directions and at the (20) and (02) sites, it was observed 
that the LDOS at the resonance energy has the strongest 
intensity directly on the Zn site, scattering from which 
is believed to be in the unitary limit. In addition, the 
LDOS at the resonance energy has nearly local minima 
at the four Cu atoms nearest-neighbor to the Zn atom, 
and has local maxima at the second- and third-nearest- 
neighbor Cu atoms. The intensity of the LDOS on the 
second-nearest-neighbor Cu atoms is larger than that on 
the third- nearest-neighbor Cu atoms. Although this un- 
expected phenomenon might indicate strong correlation 
effects, we would rather try to explain it from an alter- 
native point of view. Superconductivity in HTSC's is 
believed to originate in the Cu02 plane, and the STM/S 
tip at low bias is known to be probing such a plane clos- 
est to the surface. But the STM tip is also known to 
be separated from this plane by a BiO layer and a SrO 
layer. The BiO layer is believed to be semiconducting 
and the SrO layer, insulating. So low-energy electrons 
can go through them. However, the hard cores of the 
atoms surely can block the tunneling current from di- 
rectly going through them. We therefore postulate that 
tunneling occurs from the tip to the atoms within a small 
circular area in the Cu02 plane directly below the tip, 
except those atoms blocked in the above sense. In addi- 
tion, we postulate that the linear dimension of the small 
circular area is only a little larger than the lattice con- 
stant a, because the negative-exponential dependence of 
the tunneling current on distance implies that we need 
only include the closest set of atoms that can contribute 
to the tunneling current. (That is, they are not blocked 
by the atoms in the BiO and SrO layers closer to the sur- 
face than the Cu02 layer being probed.) Then the tip 
above a given Zn or Cu site in the Cu02 plane sees only 
its four nearest-neighbor sites in that plane. Therefore, 
we propose that the actual measured LDOS at a Zn or 
Cu site is essentially the sum of the contributions from 
its four nearest (Cu or Zn)-neighbors. The measured lo- 
cal DTC should then be related to the calculated local 
DTC by the relation: 

Gr'-(^)-EGi5-(^)- (7) 
5 

In Fig. 1^, part (a), the spatial distribution of the cal- 
culated bare [i.e., before the transformation given in 
Eq. (0)] DTC at zero bias is displaced in a three- 
dimensional plot. It includes the effects of all four 
near-zero-energy resonant [States (two per spin) localized 
around a unitary impurity.ESI The temperature is assumed 
to be at ksT = 0.02t (or T ~ O.lTc) in the calculation. 
The parameter values gi = 1.5i and /i = — 0.4i are cho- 
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sen to obtain a single near-zero-bias conductance peak. 
These values are not yet optimized, since at the present 
time we only wish to establish the essential correctness 
of our idea. In Fig. |6[ part (b), a planar "bubble plot" 
of the same data is given, where the size of each black 
dot is directly proportional to the calculated bare-DTC 
at that lattice site. (Our calculation, being tight-binding 
in nature, gives DTC values only at the lattice sites, un- 
like the observed data, which gives a continuous varia- 
tion of the DTC between the atomic sites.) In Fig. |^, 
parts (c) and (d), similar plots are given for the calcu- 
lated transformed-DTC based on Eq. (|^) . It is clear from 
comparing these plots with Fig. 3(b) and (c) of Ref.E^ 
that the calculated transformed-DTC-distribution agrees 
quite well, at least qualitatively, with the measured DTC- 
distribution at the resonant energy in that reference. 

As a crude quantitative comparison between our pre- 
diction based on Eq. (|^) and the measure data of Pan et 
al. , we have listed in Table |l| the normalized measured 
values by Pan et al. [based on Fig. 4 (a) in Ref.Ba], our 
normalized calculated bare values, and our normalized 
calculated transformed values [by using Eq. (^] , of the lo- 
cal DTC at various lattice sites near a Zn impurity [which 
is defined to be the (00) site], up to the third nearest- 
neighbor [i.e., the (20) and (02)] sites. Within each row 
of data, the normalization is such that the largest value 
becomes unity. (For the first and third rows of data, this 
occurs at the (00) site, but for the second row of data, 
this occurs at the (10) and (01) sites, because the calcu- 
lated bare value of the local DTC at the resonant energy 
vanishes at the (00) site.) We admit that this compari- 
son is only a very crude one, since our tight-binding result 
for the DTC distribution, which exists at the Zn and Cu 
sites only, and not continuously in between them, should, 
strictly speaking, be compared with some integrated re- 
sult of the measured local DTC. But in the first row of 
Table || we have only listed the measured local DTC 
values right at the Zn and near-by Cu sites, with no inte- 
gration performed. In fact, we strongly believe that the 
agreement will be much better if we do perform properly 
such an integration of the measured data before compari- 
son with our prediction, as may be seen by comparing our 
prediction given in Fig. ^, part (d), and the Fig. 3(b) of 
Ref.c3, or the three-dimensional figure on-tlie cover page 
of the March 2000 issue of Physics TodayLJ. This is par- 
ticularly so with regard to the relative heights of the local 
DTC measured at the (10) or (01) sites, and those at the 
(20) and (02) sites. The former is larger than the lat- 
ter in the data stiven in the first row of Table H, and in 
Fig. 4(a) of Ref.ta, where the data in the first row of ta- 
ble g came from. However, since at the moment we have 
not figured out the proper way to do this integration, we 
shall leave that to a future publication. (An even bet- 
ter theory should generate a continuous local DTC in a 
plane, to be directly compared with the measured data, 
with no integration needed. This possibility will also be 
looked into later.) In spite of the fact that the compar- 
ison presented in Table II is only a very crude one, we 



believe it still strongly suggests that our idea is essen- 
tially correct, although the model clearly can and should 
be improved, in order to give a local DTC distribution 
which exists continuously in a two-dimensional plane, as 
has been observed. 



B. The case of randomly distributed impurities 

In Fig. 0, we have plotted the spatially averaged DTC 
for different concentrations of impurities contained in the 
DWSC. The parameters chosen are: gi = 2t and /i = —t. 
When the super-cell contains only one impurity, because 
the spectral weight from the impurities is too small for 
this density (0.08%) of impurities, only a very small 
spatially- averaged zero-bias DTC, G{Q), appears, with- 
out a slight trace of a ZBCP. At this low concentration 
of impurities, G{y ) practically reflects the bulk d-wave 
DOS. The conductance peak outside the positive maxi- 
mum energy gap stems from the van Hove singularity. To 
increase the spectral weight of the impurity contribution 
to G(T^), one needs to increase the density of impuri- 
ties, in presumably a random distribution. But as shown 
by the dashed line in Fig. in which the randomly- 
distributed impurity density has been increased to 1.28%, 
one still does not obtain an observable ZBCP in spatially- 
averaged G{V), but only a finite spatially-averaged G{Q) 
as a local minimum. This we think is because the wave 
functions of the near-zero-energy resonant states around 
the impurities have tails along the near-nodal directions, 
the interaction between such states at different impurity 
positions (especially between those impurities not very 
far away from each other) is so large that the contribu- 
tions from these states to the DTC has spread into a 
wide band with actually lower value at the band center 
(i.e., zero bias) than at the band edges. This is con- 
sistent with the fact that in most c-axis non-STM-types 
of tunneling experiments, only a finite G(0) is observed, 
without showing a peak there of any shape and width. 

Compared with the single impurity case plotted by the 
solid line, it is clear that G(0) is enhanced at this much 
higher density of randomly-distributed impurities, show- 
ing that the spectral weight of the impurity contribution 
at this density of impurities is clearly no longer negligible. 
Yet no trace of a ZBCP is obtained. On the other hand, 
we find that at the same concentration of impurities, if 
a good portion of these impurities form [ll]-directional 
dimers, (with separation \/2a between the two impuri- 
ties in each such dimer,) but still with random orien- 
tations, a small observable ZBCP is exhibited in G{V). 
[See the short-dashed line in Fig. |^, in which 62.5% of 
the impurities form such dimers. We have also calcu- 
lated the wavefunction near the two impurities forming 
such a dimer and found that the wavefunction amplitude 
is very small along the alignment direction of the two im- 
purities, but has an oscillatory behavior perpendicular to 
the alignment direction. This highly anisotropic behav- 



6 



ior of the wavefunction may have drastically suppressed 
the interaction between the impurities. The height of 
the ZBCP depends on the total number of impurities 
forming such dimers. We also find that if some impu- 
rities form [ll]-directional trimers, the amplitude of the 
ZBCP can be enhanced even further. But we deem the 
probability for three impurities to form such a trimer in 
an actual system is very small, but the assumption that 
some impurities can form [ll]-directional dimers should 
be reasonable. Thus we propose that this is how im- 
purities can possibly contribute to some of the observed 
ZBCP's, and why often ZBCP's are not observed in c- 
axis non-localized tunneling, even in samples which have 
been deliberately introduced some substantial amount of 
substitutional non-magnetic zinc (Zn) impurities .E3 

IV. SUMMARY 

In summary, we have made detailed calculations of 
both the local and the spatially-averaged differential tun- 
neling conductance in DWSC's containing nonmagnetic 
impurities in the unitary limit. Our resjults show the fol- 
lowing: (1) Previously we have shownEj that the local 
conductance behavior near zero energy at the sites near 
a unitary impurity is sensitive to how well the particle- 
hole symmetry is satisfied, and how large is the error of 
using a semiclassical WKBJ approximation to treat the 
problem. Here we have demonstrated that conclusion 
is model independent. (2) We also find that the con- 
ductance behavior is very sensitive to the impurity con- 
figuration. For the single-impurity case, a recently ob- 
tained LDOS imaging at the resonance energy has been 
explained in terms of a model where the blocking effect 
of the atoms in the BiO and SrO layers are taken into 
account, so that the tunneling tip does not probe the 
local density of states of the Cu or Zn site directly be- 
low it, but rather essentially the sum of those of its four 
nearest neighbor sites. Furthermore, our study allows us 
to conclude that unitary impurities can contribute to an 
observable ZBCP in non-localized tunneling if and only if 
a substantial number of impurities form [ll]-directional 
dimers, (and trimers, etc.). On the other hand a sim- 
ple random distribution of unitary impurities either has 
too low a spectral weight to contribute observably to 
non-localized tunneling, or their interaction is already 
so strong that they can only produce a finite value at 
zero bias in non-localized tunneling conductance without 
giving a peak there of any shape or width. An ultimate 
test of this scenario would require: (i) the observation 
of a ZBCP in STM/S timncling in the close vicinity of 
a unitary impurity, (which is already achieved recently); 
(ii) the non-observation of a ZBCP in non-STM/S types 
of c-axis tunneling on a single-crystal sample with any 
concentration of the same type of impurities which are 
confirmed to not have formed [ll]-directional dimers, (or 
trimers, etc.); and (iii) the appearance of a small ZBCP 



in the same types of tunneling experiments when the im- 
purity distribution is established to contain such dimers, 
(or trimers, etc.). Such an ultimate test may be asking 
too much from the experimentalists, but the fact that 
this scenario is consistent with many diverse tunneling 
observations give us much confidence on it being at least 
close to the truth. 
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TABLE I. The values of d-wave order parameter and co- 
herence length for several sets of parameter values. The tem- 
perature is at T = 0.02t. 



fl/t 


91 /-t 


Ad/t 


Co/a 





2 


0.241 


1.33 


-0.4 


2 


0.229 


1.32 


-0.4 


1.5 


0.139 


2.17 


-1 


2 


0.164 


1.68 
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TABLE II. Comparison between the measured STM/S lo- 
cal differential tunneling conjdvictance at the near-zero-bias 
resonant energy by Pan et alE3 at the Zn impurity site (00), 
the nearest neighbor sites (10) and (01), the next nearest 
neighbor site (11), and the third nearest neighbor sites (20) 
and (02), (the first row), and the calculated values for the lo- 
cal differential tunneling conductance at zero bias at the same 
sites before the transformation (according to Eq. (^)) (the sec- 
ond row), and after the transformation (the third row). The 
data in each row is normalized so that the largest value is 
unity, which occurs at the (00) site in the first and the third 
rows, but at the (10) and (01) sites in the second row. 



(00) 


(10) & (01) 


(11) 


(20) & (02) 


1.00 


0.18 


0.29 


0.13 


0.000 


1.000 


0.078 


0.113 


1.000 


0.068 


0.593 


0.384 



FIG. 1. The local differential tunneling conductance as 
a function of a bias at the site directly on the site (18,18) 
(a) and (19,18) (b) in a DWSC with each supercell contain- 
ing one single impurity located at site (18,18). The impurity 
strength Uo ~ 2.5t (solid line) and lOt (dashed line). The 
other parameters /i = and gi — 2t. 



FIG. 2. The local differential tunneling conductance as a 
function of bias in a DWSC with each supercell containing 
one single impurity located at site (18,18) for fi = (solid 
line), —OAt (dashed line), and —t (short-dashed line). The 
parameter value of gi is taken to be 2t. The measure point is 
at (19,18). 



FIG. 3. The local differential tunneling conductance as a 
function of bias in a d-wave superconductor with each super- 
cell containing one single impurity (solid line) at site (18,18), 
and two (dashed line) at sites (17,19) and (18,18), and three 
impurities (short-dashed line) at sites (17,19), (18,18) and 

(19.17) aligned consecutively along the {110} direction. The 
parameters gi — 2t and /i — —t. The measure point is at site 

(19.18) . 



FIG. 4. The lowest positive eigenvalue as a function of 
the distance between two impurities aligned along the {110} 
direction in the Cu02 plane. The parameters gi — 2t and 
H = —t. 

FIG. 5. The number of near-zero energy states versus the 
number of impurities aligned consecutively along the {110} 
direction in the Cu02. The parameters gi = 2t and = —t. 

FIG. 6. Spatial distribution of the differential tunnel- 
ing conductance at zero bias: The three-dimensional dis- 
play in the whole supercell before summation over four 
nearest-neighboring Cu sites (a) and the corresponding 
two-dimensional bubble view in a small region near the im- 
purity (b); the three-dimensional display after summation (c) 
and the corresponding two-dimensional bubble view (d). The 
bubble size on each Cu site represents the conductance mag- 
nitude. The parameters g\ — 1.5t and /i = — 0.4t. 

FIG. 7. The spatially-averaged differential tunneling con- 
ductance as a function of bias in a d-wave superconductor 
with each supercell containing one single impurity [i.e., the 
impurity concentration c is about 0.08%] (solid line), a ran- 
dom distribution of impurities with c — 1.28% (dashed line), 
and 62.5% of the randomly-distributed impurities of the for- 
mer case changed to diagonal dimers (short-dashed line) . The 
parameters gi = 2t and fj. = — t. 
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FIG. 1(a) Zhu et al. 




FIG. 1(b) Zhu et al. 
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FIG. 6(c) Zhu et al. 



FIG. 6(d) Zhu et al. 
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